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We describe searches for decays to two-body charmless final states ri'ri, t/'-k'^ and jjtt" of mesons 
produced in e'^e~ annihilation. The data, collected with the BABAR detector at the Stanford 
Linear Accelerator Center, represent 232 million produced BE pairs. The results for branching 
fractions are, in units of lO"** (upper limits at 90% C.L.): B{B° -+ rj'rj) = 0.2tol ± 0.4 (< 1.7), 
B{B° riTT°) = 0.6l° j ± 0.1 (< 1.3), and B(B° ^ r/'yrO) = 0.8t° j ± 0.1 (< 2.1). The first error 
quoted is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



We present the results of searches for neutral B me- 
son decays to rj'rj, rjir^ and rj'ir'^, with a data sample 
expanded by about a factor of 2.6 over the one used 
for our previous measurements 0, 0- In the Standard 
Model (SM) the processes that contribute to these de- 
cays are described by color-suppressed tree and one-loop 
gluonic, electroweak or flavor-singlet penguin amplitudes. 
For i?'^ rj'-K^ and ^ rjiv^ the color-suppressed tree 
diagram is also suppressed by approximate cancellation 
between the amplitudes for the and for the isoscalar 
meson to contain the spectator quark, resulting from the 
mesons' isospin couplings to the quarks. Estimates of the 
branching fractions for these modes have been obtained 
from calculations based on QCD factorization 0, ^ , per- 
turbative QCD (for ^ rj^'^ir^) soft collinear ef- 
fective theory (Qj, and flavor-SU(3) symmetry The 
expectations lie in the approximate ranges 0.2-1.0 x 10^^ 
for B° rj^''^Tr°, and 0.3-2 x lO''^ for B° rj'rj. 

These decays are also of interest in constraining the ex- 
pected value of the time-dependent CP-violation asym- 
metry parameter 5/ in the decay with / = rj'K^ |7ll9lll0|. 
The leading-order SM calculation gives the equality 
S-q'K" = Sj /^k9, 1 where the latter has been precisely 
measured 'llj , and equals sin2/3 in the SM. The CP asym- 
metry in the charmless modes is sensitive to contributions 
from new physics, but also to contamination from sub- 
leading SM amplitudes. The most stringent constraint 
on such contamination in S„i iro comes from the mea- 
sured branching fractions of the three decay modes stud- 
ied in this paper 0, 0, 0| • Recently it has also been 
suggested [llllS that B° rj'n^ and B° -> 7777° can be 
used to constrain the contribution from isospin-breaking 
effects on the value of sin2a in i? ^ tt+tt^ decays. 

The results presented here are based on data col- 
lected with the BABAR detector 14] at the PEP-II asym- 
metric e+e~ collider located at the Stanford Lin- 
ear Accelerator Center. An integrated luminosity of 
211 fb^^, corresponding to 232 x 10^ BB pairs, was 
recorded at the T(4S') resonance (center-of-mass energy 
= 10.58 GeV). 

Charged particles from the e^e~ interactions are de- 
tected, and their momenta measured, by a combination 
of flve layers of double-sided silicon microstrip detectors 



TABLE I: Selection requirements on the invariant masses of 
resonances and the laboratory energies of photons from their 
decay. 

State Invariant mass (MeV) £(7) (MeV) 

7r° 120 < m(77) < 150 > 50 

ri-,-^ 490 < m(77) < 600 > 100 

7737r 520 < m(7r+7r-7r°) < 570 > 30 

77^,^^ 910 < m(7r+7r~r?) < 1000 > 100 

ij'p^ 910 < m(7r+7r"7) < 1000 > 200 

p° 510 < m(7r+7r~) < 1000 — 



and a 40-layer drift chamber, both operating in the 1.5 T 
magnetic field of a superconducting solenoid. Photons 
and electrons are identified with a CsI(Tl) electromag- 
netic calorimeter (EMC). Further charged particle iden- 
tification (PID) is provided by the average energy loss 
(dE/dx) in the tracking devices and by an internally re- 
flecting ring imaging Cherenkov detector (DIRC) cover- 
ing the central region. 

We establish the event selection criteria with the aid of 
a detailed Monte Carlo (MC) simulation of the B produc- 
tion and decay sequences, and of the detector response 
[l^. These criteria are designed to retain signal events 
with high efficiency. Applied to the data, they result 
in a sample much larger than the expected signal, but 
with well characterized backgrounds. We extract the sig- 
nal yields from this sample with a maximum likelihood 
(ML) fit. 

The -B-daughter candidates are reconstructed through 
their decays tt" 77, rj ^ iVi'y)^ V ~^ tt+tt'tt^ (ryajr), 
77' — > 77-y^7r+7r~ {rj'^^^), and additionally for rj'rj modes, 
rj' p°7 {ij'p^), where tt+tt". Tabled lists the re- 

quirements on the invariant mass of these particles' final 
states. Secondary charged pions in 77' and 77 candidates 
are rejected if classified as protons, kaons, or electrons by 
their DIRC, dE/dx, and EMC PID signatures. 

We reconstruct the B-meson candidate by combining 
the four-momenta of a pair of daughter mesons, with a 
vertex constraint if the ultimate final state includes at 
least two charged particles. Since the natural widths of 
the 77, 77', and tt" are much smaller than the resolution. 
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we also constrain their masses to nominal values m 
the fit of the B candidate. From the kinematics of T(4S') 
decay we determine the energy-substituted mass m-ES = 

yjjs — p'^g and energy difference AE = Eb — | a/s, where 
{Eb,Pb) is the B-meson 4-momentum vector, and all 
values are expressed in the T(AS) frame. The resolution 
in rriES is 3.0 MeV and in AE is 24-50 MeV, depending 
on the decay mode. We require 5.25 GeV < toes < 
5.29 GeV and \AE\ < 0.3 GcV (< 0.2 GeV for rj'ri). 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e~^e~ — *■ qq events (q = 
u,d,s,c). We reduce these with requirements on the an- 
gle 9t between the thrust axis of the B candidate in 
the T{AS) frame and that of the rest of the charged 
tracks and neutral calorimeter clusters in the event. The 
distribution is sharply peaked near | cos 0t | = 1 for qq 
jet pairs, and nearly uniform for B-meson decays. The 
requirement, which optimizes the expected signal yield 
relative to its background-dominated statistical error, is 
I cos^tI < 0.7-0.9 depending on the mode. 

In the ML fit we discriminate further against qq back- 
ground with a Fisher discriminant that combines sev- 
eral variables which characterize the energy flow in the 
event Q. It provides about one standard deviation of 
separation between B decay events and combinatorial 
background (see Fig. Qi). 

We also impose restrictions on decay angles to exclude 
the most asymmetric decays where soft-particle back- 
grounds concentrate and the acceptance changes rapidly. 
We define the decay angle 9^^^ for a meson k as the an- 
gle between the momenta of a daughter particle and the 
meson's parent, measured in the meson's rest frame. We 
require for the ry^^ decays | cos9^^^\ < 0.9 and for ry^'^Tr^ 
I cos 6'^°^ I < 0.95. For B° v'p-yVj-y the requirement is 
I '^os^docl < ^-^^ to suppress the background B K*j. 

The average number of candidates found per selected 
event is in the range 1.06 to 1.23, depending on the final 
state. We choose the candidate with the smallest value 
of a constructed from the deviations from expected 
values of one or more of the daughter resonance masses. 
From the simulation we find that this algorithm selects 
the correct-combination candidate in about two thirds 
of the events containing multiple candidates, and that it 
induces negligible bias. 

We obtain yields for each channel from a maximum 
likehhood fit with the input observables AE, toes, 
and TOi (2), the daughter invariant mass spectrum of the 
rj and/or rj' candidate. The selected sample sizes are 
given in the second column of Table Besides any sig- 
nal events they contain qq (dominant) and BB with 
h ^ c combinatorial background, and a fraction that 
we estimate from the simulation to be less than 0.2% of 
feed-across from other charmless BB modes. The latter 
events have ultimate final states different from the sig- 
nal, but with similar kinematics so that broad peaks near 



those of the signal appear in some observables, requiring 
a separate component in the probability density function 
(PDF). The likehhood function is 

N 

£ = cxp(-^y,)nE^^-^ W 

j i 3 

V,{mEslVjiAE^)Vjin'P,{^\) [-PjimD] , 

where N is the number of events in the sample, and for 
each component j, Yj is the yield of events and Vj{x^) 
the PDF for observable x in event i. For the modes 
B'^ ~^ rj'^^^rj we found no need for the BB background 
component. The factored form of the PDF indicated in 
Eq.n]is a good approximation, particularly for the combi- 
natorial qq component, since correlations among observ- 
ables measured in the data (dominantly this component) 
are small. Distortions of the fit results caused by this 
approximation are measured in simulation and included 
in the bias corrections and systematic errors discussed 
below. 

We determine the PDFs for the signal and BB back- 
ground components from fits to MC data. We calibrate 
the resolutions in AE and toes with large control samples 
of B decays to charmed final states of similar topology 
(e.g. B D{Kitit)it). For the combinatorial background 
the PDFs are determined in the fits to the data. However 
the functional forms are first deduced from fits of that 
component alone to sidebands in (toes, AE), so that we 
can validate the fit before applying it to data containing 
the signal. 

We use the following functional forms for the PDFs: 
sum of two Gaussians for ■Psig("^Es), ^sigssl^-^)' 
the sharper structures in '^^^■(toes) and Vj{mk); lin- 
ear or quadratic dependences for combinatorial compo- 
nents of Vbb qqW^k) and for Vqq{AE); and a conjunc- 
tion of two Gaussian segments below and above the peak 
with different widths, plus a broad Gaussian, for Vj{J-). 
The qq background in toes is described by the function 
— exp [—^(1 — a;^)] , with x = 2toes/-\/s and pa- 
rameter ^. These are discussed in more detail in and 
some of them are illustrated in Fig. 2] 

We allow the parameters most important for the de- 
termination of the combinatorial background PDFs to 
vary in the fit, along with the yields for all components. 
Specifically, the free background parameters are most or 
all of the following, depending on the decay mode: ^ 
for Toes, linear and quadratic coefficients for AE, area 
and slope of the combinatorial component for m^, and 
the mean, width, and width difference parameters for J^. 
Results for the yields are presented in the third column 
of Table HTl for each sample. 

We test and calibrate the fitting procedure by apply- 
ing it to ensembles of simulated qq experiments drawn 
from the PDF into which we have embedded the expected 
number of signal and BB background events randomly 
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TABLE II: Number of events A'' in the sample, fitted signal yield Ys in events (ev.), measured bias, detection efficiency e, 
daughter branching fraction product (]~[ Bi), and measured branching fraction B with statistical error for each decay chain, and 
for the combined measurements the significance S (with systematic uncertainties included), branching fraction with statistical 
and systematic error, and in parentheses the 90% C.L. upper limits. The number of produced BB pairs is (231.8 ± 2.6) x 10®. 

Mode iV (ev.) Ys (ev.) Bias (ev.) e (%) UB, (%) S (a) B (lO"'^) 



^777r7r^37r 


539 


2 n+^'i 

^•^-2.0 


1.9 ±1.0 


13.8 


3.95 




1+2'' 








1448 


9 1+3.5 
^•-'--2.2 


0.7 ±0.4 


22.3 


6.89 




4+'° 








8268 


-8.61?- 

1 c + lO.S 


0.0 ±0.4 


14.9 


6.67 




Q Q + 3.8 
-•3-8-3.0 






Vp-yVjt 


16861 


0.0 ±0.5 


21.8 


11.63 




0.2lli 






v'v 












0.4 




± 0.4 


«l-7) 




2334 




1.2 ±0.7 


16.3 


22.6 




1 1+1.0 

^■^-o.s 








5493 


6.51^^6^ 


1.2 ± 0.8 


20.7 


39.4 






















1.3 




± 0.1 


«1.3) 




3663 


7.9+_ti 


1.2 ±0.6 


17.5 


17.5 


1.4 




± 0.1 


«2.1) 



extracted from the fully simulated MC samples. We find 
biases of 0-2 events, somewhat dependent on the signal 
size. The bias values obtained for simulations that repro- 
duce the yields found in the data are given in the fourth 
column of Table m 

In Fig. n] we show, as representative of the sev- 
eral fits, the projections of the PDF and data for the 
r/^^^7r" sample. The goodness-of-fit is further 
demonstrated by the distribution of the likelihood ratio 
C'bkg] for data and for simulation gener- 
ated from the PDF model, shown for the same decay 
mode in Fig. [21 We see good agreement between the 
model and the data. By construction the background is 
concentrated near zero, while any signal would appear in 
a peak near one. 

We determine the reconstruction efficiencies, given in 
Table ^1 as the ratio of reconstructed and accepted 
events in simulation to the number generated. We com- 
pute the branching fraction for each channel by subtract- 
ing the fit bias from the measured yield, and dividing the 
result _by the efficiency and the number of produced BB 
pairs We assume equal decay rates of the T{AS) to 
B^B^ and B^B^ . Table ITU gives the numbers pertinent 
to these computations. The statistical error on the signal 
yield or branching fraction is taken as the change in the 
central value when the quantity — 2 ln£ increases by one 
unit from its minimum value. 

We combine results where we have multi- 
ple decay channels by adding the functions 
-21n{[£(S)//;(So)] <^G(S;0,ct')}, where 6o is the 
central value from the fit, a' is the systematic uncer- 
tainty, and ®G denotes convolution with a Gaussian 
function. We give the resulting final branching fractions 
for each mode in Table ITU with the significance, taken as 
the square root of the difference between the value of 
— 2 ln£ (with additive systematic uncertainties included) 
for zero signal and the value at its minimum. The 90% 
C.L. upper limits are taken to be the branching fraction 
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FIG. 1: Plots of the — > ry^^^vr'' data distribution projected 
on each of the fit variables: (a) AiJ, (b) tties, (c) t)' mass, and 
(d) T. The solid line represents the result of the fit, and the 
dashed line the background contribution. (The absence of 
signal here nearly hides the dashed curve.) The dotted line il- 
lustrates the expected shape for signal, normalized arbitrarily 
to the data. 



below which lies 90% of the total of the likelihood 
integral in the positive branching fraction region. 

The systematic uncertainties on the branching frac- 
tions arising from lack of knowledge of the PDFs have 
been included in part in the statistical error since most 
background parameters are free in the fit. For the signal, 
the uncertainties in PDF parameters are estimated from 
the consistency of fits to MC and data in control modes. 
Varying the signal-PDF parameters within these errors, 
we estimate yield uncertainties of 0-2 events, depend- 
ing on the mode. The uncertainty from fit bias (Table 
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Ty^^^vr". The open circles represent an arbitrarily large sim- 
ulated signal component, the solid points represent the data, 
the solid histograms are from toy samples of background 
(shaded) and background plus signal (white, barely visible 
in the rightmost bins, given the small signal yield). 



value of the CP asymmetry Sf'va. relation to sin2/3 for the 
decay i?° -q'K^ llfll ■ Using the method proposed 
by Gronau et al. jlflj, we estimate that our results will 
provide approximately 20% improvement of the predic- 
tion for the contribution of the color suppressed tree am- 
plitude in — > Tj'K^ decays. This translates into a 20% 
reduction of this theoretical uncertainty in Sj^ij^o. We 
find a similar improvement in the correspondin g un cer- 
tainty of sin2a measured with _B ^ tt+tt^ decays |l2Lll3j . 
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chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the computing 
organizations that support BABAR. The collaborating 
institutions wish to thank SLAG for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DEC (Germany), 
INEN (Italy), EDM (The Netherlands), NER (Norway), 
MIST (Russia), and PPARC (United Kingdom). Indi- 
viduals have received support from CONACyT (Mex- 
ico), Marie Curie EIF (European Union), the A. P. Sloan 
Foundation, the Research Corporation, and the Alexan- 
der von Humboldt Foundation. 



Inj) includes its statistical uncertainty from the simulated 
experiments, and half of the correction itself, added in 
quadrature. Similarly we estimate the uncertainty from 
modeling the BB backgrounds by taking half of the con- 
tribution of that component to the fitted signal yield, 
0.2-1.2 events. These additive systematic errors are dom- 
inant for these modes with little or no signal yield. 

Uncertainties in our knowledge of the efficiency, found 
from auxiliary studies, include 0.8% x Nt and 1.5% x N^, 
where Nt and N-y are the number of tracks and photons, 
respectively, in the B candidate. The uncertainty in the 
total number of BB pairs in the data sample is 1.1%. 
Published data provide the uncertainties in the B- 
daughter product branching fractions (0.7-3.9%). The 
uncertainties in the efficiency from the event selection 
are about 1% . 

After combining the measurements we obtain the cen- 
tral values and 90% C.L. upper limits for the branching 
fractions: 

B{B° ?/r/) = (0.2t°-^ ± 0.4) x 10"^ (< 1.7 x 10"^), 
B{B^ riTT°) = (0.6i[j:^ ± 0.1) X 10"^ (< 1.3 x lO^^), 
and 

B{B^ r]'TT°) = {O.Stoi ± 0.1) X lO-^ (< 2.1 x 10-^). 

We find no evidence for these decays, and our upper lim- 
its represent two to three-fold improvement over the pre- 
vious measurements P, 0, ESlJ ■ The range of sensitivity 
of these measurements is comparable to the range of the 
theoretical estimates. 

These results can be used to constrain the expected 
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